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Using flow equations, we derive an effective quark-quark interaction and obtain 
the coupled set of gap equations for the condensates of the CFL phase of massless 
Nt = 3 dense QCD. We find two different sources of the infrared cutoff in magnetic 
interaction. When the penetration depth of the magnetic field inside the supercon- 
ductor is less than the coherence length of quark-quark bound state, our results 
for the gap agree with those of Son 4 . In the other case, we obtain parametric 
enhancement of the gap on the coupling constant. 



1. Introduction 

Quark matter at sufficiently high density is a color superconductor. It has 
been known for some time 1 , however the value of the gap was too small to 
be of any practical interest. Recent revision of the subject brought much 
larger gap 2 than previously thought. The consistent numerical value of the 
gap has been obtained using phenomenological four-fermion interactions 3 
and from QCD one-gluon exchange 4 , though with different parametric 
dependence on the coupling. Many microscopic calculations of the gap in 
Nf = 2 5,6 and in Nt = 3 7 have been done after that. 

Microscopic studies rely on the fact that due to asymptotic freedom 
QCD becomes a weakly interacting theory at high densities and perturba- 
tion theory can be used. However, high density quark matter cannot be 
described as a simple Fermi gas. As known from the BSC theory of super- 
conductivity, the Fermi surface is unstable to an arbitrarily weak attractive 
interaction. In QCD attractive interaction is mediated by one-gluon ex- 
change between quarks in color antitriplet state. 
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In the language of renormalization group, BCS instability with respect 
to Cooper pairing is associated with breaking down the conventional per- 
turbation theory in the vicinity of the Fermi surface. This signals formation 
of a gap which regulates the infrared divergent behavior. It is in full anal- 
ogy with electron superconductivity in metals. The difference between the 
two is that the attractive interaction via gluon exchange is already present 
in QCD rather than being induced by the solid state lattice as in metal. 
Therefore, properties of a gluon propagating in high density quark mat- 
ter are essential in calculating the color superconducting gap. It has been 
suggested by Son (the first reference in 4 ), that the gap is dominated by 
magnetic long range exchanges, and that the infrared behavior is regu- 
lated by dynamic screening in quark plasma. We show that there are two 
characteristic scales in superconducting medium: the penetration depth of 
magnetic field inside the superconductor, 5, and the coherence length of 
the bound quarks (or size of the qq Cooper pair), £. The prevailing en- 
ergy 1/5 or l/£ regulates the long range behavior of magnetic field in color 
superconductor. 

The purpose of the present work is to obtain an effective microscopic 
low energy theory of quasiparticles and holes in the vicinity of the Fermi 
surface. We shall argue that depending on the state of superconductor, 
there are two different sources of the infrared cutoff in magnetic interaction. 
At small temperature, the superconductor is of a Pipperd type, S <C £, 
and magnetic gluon is damped by the quark medium according to Landau 
damping mechanism. As we increase the temperature, the quark plasma 
'dilutes' and the penetration depth of magnetic field grows. It turns out that 
usually the superconducting condensate melts slower than the penetration 
depth grows 12 , corresponding to London type of superconductor, £ <C 6. 
In this case magnetic gluon scatters over diquark Cooper pairs and damped 
by the superconducting medium. Magnetic interaction is regulated by the 
inverse coherence length, or roughly by the superconducting gap. In these 
two limiting cases we obtain for the gap the same parametric dependence 
on coupling constant, though different numerical factors. 

We apply flow equations to the Coulomb gauge QCD Hamiltonian with 
Nf = 3 at nonzero quark density. Flow equation method 10 , which is a 
synthesis of the pcrturbativc renormalization group and the many-body 
technique, has been successfully applied to the problem of electron super- 
conductivity in metals 11 . The main idea of the method is to eliminate 
interactions which mix states with different particle content (or number of 
particles) and obtain an effective theory which conserves particle number 
in each (Fock) sector. This procedure is performed in sequence for matrix 
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elements with large energy differences down to more degenerate states. In 
high density QCD, all nonabelian contributions are suppressed, therefore, 
the quark-gluon coupling is the only term which mediates hoping between 
states of high- (away from Fermi surface) and low- (near Fermi surface) 
energies. Eliminating quark-gluon coupling by flow equations we virtually 
move towards the Fermi surface, in analogy to perturbative renormalization 
group scaling. Decoupling different energy scales corresponds to integrat- 
ing out fast moving modes - quarks in the Dirac sea, and obtaining an 
effective theory for slow modes - quarks close to Fermi surface. As men- 
tioned above, in this scaling down we encounter two characteristic scales: 
penetration depth of magnetic field and the coherence length of qq bound 
state. Depending which scale is bigger, we obtain different mechanisms of 
the infrared cutoff in magnetic interaction and corresponding gaps. 



2. Effective microscopic Hamiltonian for color 
superconductivity 

We start with the Coulomb gauge QCD Hamiltonian,V-A = 0, with Nf = 3 
at nonzero quark densities given by 

H = H + Hi = H + H mst + H dyn , (1) 

where H is the free Hamiltonian, H inst is the instantaneous interaction 
describing static properties, and Hdyn is the dynamical interaction involving 
the gluon propagation. The free part is the kinetic energy of quarks and 
gluons 

H a = J dxtp(x) {—i-f ■ V — ^70 + m) ip(x) 

+ TrJ dx(U 2 (x)+B 2 A {x)) , (2) 

where the non-abelian magnetic field is B = Bi = VjAk — VkAj+g[Aj, A^], 
and its abelian part is represented by Ba- The degrees of freedom are the 
transverse perturbative gluon field A = A a T a {A = A PT ), its conjugate 
momentum Et, and the quark field in the Coulomb gauge. The instanta- 
neous interaction is given by 

H inst = -\j dxdy^(x) l0 T a ^(x)V mst (\x - y\Myh T a yj(y) , (3) 

where the leading order kernel is a Coulomb potential defined, together 
with its Fourier transform, by 

2 

V ins t(r) = -— , Wg) = -Hj, ( 4 ) 
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with a s = g 2 /4:n. The dynamical interaction includes the minimal quark- 
gluon coupling, V qg , and the non-abelian three- and four-gluon interactions, 
V gg , i.e. H d yn = Vqg + Vgg, where 

Vqg = -g J dxi){x)~f ■ A(x)i>(x) 

V gg = Trj dx (B 2 (x) - B\{x)) . (5) 

As mentioned in the introduction, at high quark densities diagrams in- 
cluding the non-abelian interactions are suppressed, henceforth V gg is not 
considered. 

Following the work of Alford, Rajagopal and Wilczek, 3 , we assume, 
due to quark interactions, the diquark condensation in scalar < ip T C^tp > 
and pseudoscalar < tp T Cip > channels, or equivalcntly < bb >^ and 
< ww >t^ 0. In other words, pcrturbative vacuum |0) is not a ground 
state for the system, instead it is the BCS vacuum \Q) containing conden- 
sates of diquarks. The Fock space is constructed from this vacuum using 
quasiparicle operators and d) 

#*0 = E / j^Hk,s)b(k,s) + v(-k,s)dH-k,s)}e ikx 

n(*) = -<E / v^s\I^Hk,a)-aH-k,a)]e ikx , (6) 

where 6|f2) = d|0) = 0, and the gluon part has trivial vacuum o|0) = 
with u>(k) — \k\. All descrete numbers (helicity, color, and flavor for the 
quarks and color for the gluons) are collectively denoted as s and a, re- 
spectively The gluon operators a = a IA (k) = J2x=i 2 X)a A (k, A) 
are transverse, i.e. k ■ a A (k) = and the polarization sum is Djj(k) = 
5^ A=1 2 e/(fc, A)ej(fc, A) = Sjj — k]kj. In the massless basis, m = 
0, u{k) = v(k) = (L(k),R{k)) with the left Weyl spinor L(k) = 
(- sin(6(fe)/2) exp(-^(fe)), cos(6(fe)/2)). 

Then, in the mean field approximation, the Hamiltonian Eq. (1) can be 
written in terms of creation/annihilation operators as 

/dk f dk 

^ (k - m) bHky a b(ky a + J — {p - k ) bHky a b(ky a 

/dk ( dk 

(k + n) dHky a d(ky a + J — k at(fe)a(fe) 
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+ f A A^(p)T^-W) (j(p^(-p)j + ^)^(-p)j) + <■■<■. 
W(k,p)e l 



(2tt) 3 

— -wrfc „^<*(fc) fi -**(P) 

(2tt) 3 (2tt) 3 4 
x 6t (fe )i KT ^ 6(p )^t ( _ fe) i T A 6( _ p) ^ + . . . 

x tt( fc) ir^(p)' fl ^ fc _~^ ( u t(fc)^(p)) + . . . , (7) 

we explicitly display color (a, /3), flavor (i, j) and polarization / indices. The 
first three terms describe kinetic energies of particles, holes and antipati- 
cles, respectively. The diquark condensate A and the effective quark-quark 
interaction W are unknown parameters, and dots show that there are other 
possible terms which we ignored. 

Color-superconducting condensate A^ 7 is a N c x N c matrix in funda- 
mental color space (a, 7 = 1, N c ), and a Nf x Nf matrix in flavor space 
= 1, Nf). It can be parametrized for Nf = 3 as 3 

3 



A 8 (p) + ^Ai(p) 



^ + iA 1 (p)On , (8) 



where Ai and A 8 are the gaps in the singlet and octet channels, respectively. 
We will also use the parametrization in terms of symmetric, A( 6 6 ) , and anti- 
symmetric, A(3 3), in color and flavor gap funtions, A^ 7 = A( 6j6 )(<^# 7 + 
$y<%) + A(3,3) (^a^7 — The connection with the eigenvalue gaps is given 
by A ( 3 i3) = 1/2(A 8 - l/4Ai), and A (6;6) = 1/2(A 8 + l/2Ai). 

Replacing color a and flavor i indices with a single color-flavor index p, 
A£? 7 is diagonal in the CFL basis 9 

b(ky a = E^ 5 ( fc ) p ' ( 9 ) 

where A p are the Gell-Mann matrices for p = 1, ...,8 and A 9 Q = a/2/3<5^ 
for p = 9, and strictly speaking we should have used different letter no- 
tation for the CFL operators b p . Then, A^{p)b(p)^b(-p) 3 s T^T^ s = 
f /2 V A p (p)b(p) p b(-p) p , and the Hamiltonian Eq. (7) is given in the 
CFL basis by 

H = Y,\k-fi\ b\{k)b p {k) + ]T(fc + n) 4(fc)d p (fc) + ^ fc ot(fc)o(fc) 
fe.p fe,p k 

+ \ £ A p (p)e-^(P) (& P (P)M-P) + 4(P)4("P)) + c - c - 
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+ 3 £ {k,p)±K,^>&>g s T%TfrW k h-«W 
k,p,p, P ' 

x bl(k)b p ,(pp p (-k)b p ,(-p) + ... 
k,p,p, P ' 

x bl(k)b pl (p)^-^(uHk)a I u( P )) + ... . (10) 

where we used notation J2k = / dk/(2ir) 3 . As in 3 , we change basis to 
creation/annihilation operators y and z for quasiparticlcs (quasiholes) and 
quasiantiparticles, respectively, 

y p (k) = cos(e»(fc))^(fc) +sin(e»(fc))exp(i^(fc))6+(-fc) 

z p (fc) = cos(e^(fc))6"(fc) + sin(0*(fc)) exp(i^(fc))6+(-fc) , (11) 

with j/ p (fe)|0) = z p (k)\0) = 0. In order to absorb the condensate term 
into a new free Hamiltonian, we choose 3 cos(29^(fc)) = \k - fj,\/E { p \k), 
cos(26*(fc)) = \k + n\/E { p + \k), and sin(2e»(fc)) = A p (k)/ £ p _) (k), 
sin(29*(fc)) = A p (fe)/4 +) (fc); $(k) = cj>(k) + it, £*(fc) = -#fc). As 
the result of this transformation, a new free Hamiltonian is given by 

Ho = bl(k)b p (k) + 5>< + >(fc) 4(k)d p (k) 

k,p k,p 

+ ^w(fc) at(fc)a(fc), (12) 

where £ p _) = V( fc - m) 2 + A p (fc) 2 , 4 +) = V( fc + A*) 2 + A p (fc) 2 , and 
w(k) = fc. The effective Hamiltonian Eq. (10) is given by 

H = H + V qq + V qg , (13) 

with transformations of Eq. (11) made in the diquark interaction V qq and 
the quark-gluon coupling V qg , which are the final two terms in Eq. (10). In 
the next section, our aim is to determine the unknown parameters of the 
effective Hamiltonian Eq. (13) using flow equations. 



3. Flow equations 

Flow equations are written for the unknown functions W pp , g pp , A p of 
the effective Hamiltonian Eq. (13). Calculations are performed using y, z 
variables and the vacuum |0). Similar calculations have been done in 8 . 
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The first order flow equations dV qg /dl = [t],Ho], i] = [Ho,V gg ] eliminate 
the quark-gluon coupling V qg . The generator of the transformation is given 

by 

r,= ]T V pp '(k,p,k-p)^\ P a T^ 
k,p,p, P ' 

x bt(fe) V (p) a ^ fc _-P) ( M t (fc)a J M(p)) ; (14) 

then for the generator rf p and coupling g pp functions flow equations are 
written as 

= -(E p (k) E p ,{p) - P) f 9 PP 

V pp ' = (E p (k) - E p ,(p) - u(k - p)) g p p' . (15) 

The solution of Eq. (15) reads 

gp"' {I) =g(0)e X p(-(E p (k)-E p ,(p)-u(k-p)f l) , (16) 

where g(0) = g is the bare coupling constant. As I — > oo the coupling is 
eliminated as long as the states in the exponent are not degenerate. 

Using the generator Eq. (14) and the quark-gluon coupling Eq. (10), 
we write the second order flow equations for the effective quark-quark 
interaction, dV qq /dl = [77, V qg ]t wo -body, and for the diquark self-energy, 
dYi q /dl = [r/, Vqg]one-body Flow equation for the diquark interaction is 
given by 

^p = - ^(fc)6t(-fc) 6 ^(p)V(-p)^ a A^iA?;A^T^ 

k,p,prho' 

(u^k)a\(p))(u^-k)a J u(-p)) D 2 I J^-^ (17) 

{rf 1 (k, p,k- p)gp' p {p, k,k-p) + V p ' p (p, k,k- p) 9 pp' (k, p,k- p)) . 

Integrating the flow equation Eq. (17) over I = [0, 00), with the final value 
V qq — V qq (l — 00), and adding the instantaneous interaction Eq. (3), we 
get 

k,p,p, P ' 

(v p p'(k,p)(uHk)a I u(p))(u\-k)a J u(-p))D IJ (k - p) 
+V{k,p){u\k)u{p)){u\-k)u{-p))) , (18) 
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with 
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V" (k,p) = - 



2 (EP{k)-EP'(p)) 2 + oj M (k-p) 2 



V(^p) = ^r f, 1 , 2 , (19) 
2 We(k - pY 

where we have used solutions for rf p and g pp , Eqs. (15,16). Energies for 
magnetic and electric gluons, lom and loe respectively, are specified further. 
Keeping only left-left components, matrix elements are given by 

{u\k)a I u(p))(u\-k)a J u{-p))D IJ (k - p) = 
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1 j ^ 2 + 2 (k-p) 2 J 

{u\k)u{p)){u\-k)u{-p)) = (-lJe-We^l (^^) - ( 20 ) 

where we used Fierz transform (u(fc)O j4 u(p))(fl(— k)0 B u{— p)) = 
eEcDfcS^c(p)O c u(-p))(u(-k)0 D u c (k)), where e = -1 for 
O a ,O s = 7 M , and Cw T = uc, u T C = uc- For the above matrix ele- 
ments the coefficients / in Fierz transform O a x O b — > O c x O 15 are given 
by 70x70 — > 1/4(1 x I-75 x 75) + 1/4(7^ xct^", and7 4 X7 J -» <r?l/4(lxl- 
75 x 75) + ^1/40-^ x o pv + 1/4(0^ x oJ> + t7„i x <7« - x <7« - a,,, x <^>), 
where /x, f = (0, i = 1,2,3), 0^ = 1/2[7 M ,7„]. The diquark condensates 
are 1 x 1-pseudoscalar (ucu), 75 x 75-scalar (ucju), and x er^-vector 
{uc^iivu). Performing algebra with A-matrices in Eq. (18), we get 

^ = -E(- 1 ) c ^ (p)e40(fe) 

k.p 

(4(fc)4(-fc)fe 1 (p)fe 1 (-p)f 81 (fc,p) + 6t( fe)fo t ( _ fc)fe8(p)68( _ p) ^i8 (A , jp) 

- bl(k)bl(-k)b 8 (p)b s (-p)2F ss (k 7 p)) , (21) 
where we introduced 

F PP ' = V&(k,p) {^^j + V inst (k,p) (^P^j , (22) 

and 

v pp' (k v) = 2 jL I 

v dyn\ >P) 3 (EP(k)-EP'{p)) 2 +LO M (k-p) 2 

V mst (k,p) = \- 1 . (23) 
3 uj E {k-p) 2 
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Dynamically generated by flow equations interaction describes magnetic 
gluon exchange, and the instantaneous interaction -electric gluon exchange. 
Comparing Eq. (21) with Eq. (10) for V qq , the effective quark-quark inter- 
action is given by 

W»»' (k,p) = -F pp ' (fe, p) = -V d p y i(k, p) - Vf n i (k, p) , (24) 

in the collincar limit k ■ p — 1. This interaction is attractive in the singlet 
and octet channels. Note that W 11 = 0, and condensation in singlet channel 
is driven by W 81 . 

Dynamical magnetic interaction, Eq. (23), has the form — l/(q 2 +5E 2 ) 
instead of — l/(q 2 — SE 2 ), that is produced by the equal time perturbation 
theory, where SE is the enegry difference of in- and out-going quarks. The 
latter has the pole, therefore usually SE is neglected near the Fermi surface. 
Our interaction is regular, and SE will play an important role. 

In order to incorporate effects of the dense quark medium, we include 
polarization operators for electric and magnetic gluons, modifying the gluon 
single energy as done in the work of Pisarski and Rischke 5 

1 If (k-p) 4 



u M (k -p) 2 2 V (fe - pf + M 4 (E(k) + E(p)) 2 ) 

(k-p) 4 



+ 



(k - pf + M A (E(k) - E{p)) 2 ) 
1 



UE{k — p) 2 (k—p) 2 + 3m 2 



(25) 



a 

where M 2 — (37r/4)m 2 , m 2 = N / g 2 fx 2 / '{&ir) 2 . Eq. (25) is used in the 
complete diquark interaction Eq. (23). 

Quark self-energy has b^b and bb + tftf components associated with 
normal and anomalous propagation, respectively. For the anomalous prop- 
agation self-energy the flow equation is given by 

d^qq _ _ 1 X( o x p 1 xp' x p' rpA rpA 

k,p,p,p> 

(< &t(fc) 6 t(_ fc) > v(p)M-p)+ < V(p)M-p) > b\(k)bl(-k)) 

{u\k)aU{p)){u\-k)a J u{-p)) D ^-^ 

(if' (fe, p, k - p)g p ' p {p, k,k-p) + V p ' p (p, fe, fe - p)g pp ' (fe, p, k - p)) 

+ (Terms ~< dd >, < >) . (26) 

We neglect the contribution of antiparticles, since it is regular near the 
Fermi surface, < dd >^ sin(29 z ) <~ A/2/i, and does not change much the 
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condensate of particles. The anomalous propagator for particles is given by 

< &t(fc)6t,(-fc) > = -I sin (2e,(fe))ex P (-^(fe)), (27) 

where the average is calculated in |0) vacuum. 

Integrating over I = [0, 1/A 2 ], where A is the UV cut-off, and adding the 
self-energy instantaneous term comming from normal-ordering Eq. (3) in 
the |0) vacuum, we get 

(y ppl (k,p)(u\k)a I u(p))(u^(-k)a J u(-p))D IJ (k - p) 
+F(fc,p)(u t (fc)u(p))(u t (-fc)u(-p))) R pp '(k,p;X)+c.c, (28) 

where V p ' p and V are given in Eq. (19), and the the UV regulator is given 

by 

R pp '(k,p;\) = C ^{-[{E p {k)-E p '{p)f+uj M {k-p) 2 ]/\ 2 ) , (29) 

Change in sign appear due to adopting conventional scaling up as in the 
perturbation theory, instead of scaling down in flow equations. Performing 
calculations with A-matrices, we have 

Z gq = £(-l)e-^*V<« fe ) (&i(p)6i(-P) < bl(k)bl(-k) > F 81 (fc,p) 
k,p 

+MPM-P)[< b\(k)b\(-k) > F ls (k,p)- < bl(k)bl(-k) > 2F ss (k,p)}) 
x R pp '(k,p;\) +c.c. . (30) 

Comparing Eq. (30) with Eq. (10) for the gap functions, we obtain the 
system of gap equations 

A 1 (p) = 8G 81 (p) 

A 8 ( P ) = G 18 (p)-2G 88 (p), (31) 

where 

, i r flu i , 

G pp (p) =--J —-sm{2Q p {k))FP>> (k,p)R pp (k,p;X), (32) 

with F pp is given by Eq. (22), and R pp by Eq. (29). 
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4. Solving the system of gap equations 

We solve approximately the system of gap equations Eq. (31) in two lim- 
iting cases of the Pipperd and London type superconductor. Magnetic 
interaction Eq. (23) is regulated by the term generated by flow equations, 
with dispersion q ~ E, and by the magnetic gluon polarization operator, 
Eq. (25), with dispersion q ~ E 1 /^ . When the first term dominates, super- 
conductor is of London type, the second term-Pipperd type. 

We take instead of the smooth regulator R, Eq. (29), a sharp cut-off, 
\k — /i| < <5. Near the Fermi surface the integration measure is given by 
J dk = 2nfi 2 J_ s d(k — fi) J_ x dcos9, and since the integral is even in (fc — fi) 
one has Jg — » 2 J"*. In the denominator of F, Eq. (22), we approximate 
(fc-p) 2 = 2 M 2 (1 -cos 6»). 



Performing the ^-integration, we get for the Pipperd superconductor: 
where 

S = ; ^ = 32?r (j^) ' 5M+£ = 256?r4 (^?) 5/2 ' (34) 

and superscript M and M + S means that only magnetic or magnetic and 
electric components are taken. For the London superconductor: 

lj f S d(k-,,) A ,,,1, / i>V 



^w-i^ w^ k w)-w '' (35) 



where 



1/2 

2 , b M+E = 4tt ( — ] . (36) 



Substituting G pp Eqs. (33,35) into Eq. (31), we have the system of Eliash- 
berg type of equations. Different effective coupling g is obtained because 
of different dispersion law in the two cases. The coupling g defines the 
exponent of the solution for the gap A and b the preexponential factor. 

If we are interested only in the correct exponential dependence, we can 
neglect p dependence by the quark energies in the anomalous quark propa- 
gator (sin factor in G) and in the gluon propagator (under the In). Indeed, 
when estimating the momentum-independent gap, double logarithm arises 
from the region A <C k <C \i where the integral reduces to 
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with the solution A ~ 6^exp(— \/2/g); i.e. double logarithm responsible 
for the exponential solution does not depend on the quark energies. To 
calculate the preexponential factor correctly, one should make the rescaling 
in the integral as suggested by Schafer 7 . 

In order to convert the integral gap equation into differential one, we 
split the logarithm in G pp as suggested by Son 4 

ln (^))^-^ + ln (^))^" P) ' (38) 
that gives the same expression in two cases, Eq. (33) and Eq. (35). Coupled 
gap equations Eq. (31) decouple for the antitriplet and sixtet in color and 
flavor gaps, introduced after Eq. (8), 



Introducing the variable, as by Pisarski and Rischke 5 , 

x = ln ( z 2b ^r^ ) > ( 4 0) 



2 



/ 6/i \ f (p - p) d{k - n) 



the integratiom measure is given by d(k — fi)/E(k) = dx/x, and the integral 
equations, Eq. (39), reduce to differential equations 

^A(3 i 3 ) +5 A (3i3) =0 
dp 1 -2 

— A (6i6) - |-A (6i6) =0, (41) 

with initial conditions at the Fermi surface dA/dx(x = x ) = and A(x = 
xo) = A ; also away from the Fermi surface A(x = 0) = 0. The negative 
sign in the second equation means the repulsion in the sixtet channel. As 
function of momentum there is a trivial solution in the sixtet channel. In 
the antitriplet channel, the solution reads 

A (3,3) = A sin(gx) , A = 26/icxp(-^) , (42) 
where g and b are given in Eqs. (34,36). 
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5. Conclusions 

Using flow equations, we derived an effective quark-quark interaction and 
obtained the coupled set of gap equations for the condensates of the CFL 
phase of massless Nf = 3 dense QCD. Diquark interaction, generated dy- 
namically by flow equations, is a long-range magnetic gluon exchange reg- 
ulated by two different sources in the infrared region. One term describes 
the retardation effects of a magnetic gluon caused by Landau damping in 
dense quark-gluon plasma (normal phase). At small temperature, far be- 
low the melting of superconducting condensate T C To, Landau damping 
is the dominant mechanism in the infrared, that corresponds to the Pip- 
perd type of superconductor. The other term describes retardation due to 
propagating of a magnetic gluon in a superconducting matter, i.e. through 
the multiple scattering of a gluon at diquark Cooper pairs. This mechanism 
is dominant at temperatures close to a melting point of superconducting 
condensate, T — T <C T , and corresponds to the London type of super- 
conductor. 

We obtain approximate analytical solutions of the gap equations in these 
two limiting cases. The dominant contribution to the condensate comes in 
the color antitriplet, flavor antitriplet channel. The color sextet, flavor sex- 
tet contribution is small but non-zero. In the color and flavor antisymmetric 
channel, we obtain parametric enhacemcnt of the London type condensate 



The same conclusion holds for Nf = 2. Numerically, the (3, 3) codensates 
are almost the same for the London and Pipperd type 2SC superconductors, 
Nf = 2. However, the condensate of the London type superconductor is 
sufficiently larger than of the Pipperd superconductor in the CFL phase, 
Nf = 3. This has important implications, which regime is energetically 
favorable and will actually be realized in neutron stars. 

Independence of the gap on number of flavors for magnetic field and slow 
dependence on Nf for the sum of magnetic and electric fields in London 
type superconductors works in favor that the CFL phase will be observed 
in neutron stars, and that there is no window for the 2SC phase 13 . 

Our calculations have been done at zero temperature. Rigorous cal- 
culations of the condensates which include calculations of the penetration 




(43) 



versus the Pipperd type condensate 



A ~ exp( 




(44) 
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depth of magnetic field and coherence length at non-zero temperature are 
necessary. Qualitatively it will not affect our results, since our conclusions 
rely on similar trends known in metal superconductors, which are derived 
from Landau-Ginzburg theory and are proven experimentally 12 . 

We included all possible condensates, < qCq >, < qC-f^q > and < 
qCa^q >, which do not mix left and right components, that accounted 
for the angle dependence in the gap equations. It is interesting to add the 
instanton-induced 't Hooft interaction and to consider simultaneously the 
chiral condensation in qq channel. Generalization to a non-zero strange 
quark mass, and adding kaon condensates is important for neutron star 
phenomenology 
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